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INTRODUCTION 


THE distribution of plant species in plant communities has been studied 
from several different points of view, but few mathematical analyses of 
the problem have been made. The labor in the complete mapping of all 
the individuals of a single species in a community usually precludes this 
type of presentation and less satisfactory schemes must be used. 

The usual method of obtaining quadrat frequency data on a single 
species is the sampling of a representative area of the community by a 
number of scattered areas (quadrats) and the determination of the num- 
ber of plants of that species in each quadrat. If the plants are distri- 
buted completely at random over the area and they have no influence on 
each other it can be shown statistically that the distribution will be of 
the Poisson type (see below). This was established for a number of 
species in certain grassland communities by Blackman (1935) and by 
Ashby (1935). The presence of clones, rhizomes, and runners and other 
clumping factors, however, vitiates the fundamental assumptions and 
under these conditions the data do not follow a simple Poisson law. 
Among the workers who have studied the general problem of nonrandom- 
ness are Rommel (1930), McGinnies (1934), Ashby (1935, 1948), Blackman 
(1935, 1942), Clapham (1936), Fracker and Brischle (1944), Cole (1946), 
Archibald (1948), Whitford (1949), Thomas (1949), and Curtis and Mc- 
Intosh (1950). Many of these investigators proposed measures of non- 
randomness (see below), but the general problem of obtaining a measure 
of the size of the clumps was not taken up by them. 
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Two similar lines of approach to these nonrandom distributions are 
afforded by the use of "contagious" distribution functions such as those 
of Neyman (1939) and Thomas (1949), where the mathematical model 
assumes that primary centers ("clusters") are distributed at random 
and that individual plants are then distributed, again at random, about 


these primary centers. These models have proved more or less success- 


ful in representing certain plant-distribution data (Archibald, 1948; 
Thomas, 1949). A typical example of such an excellent agreement is 
provided by these authors for the distribution of Armeria maritima in 
an English marsh (Table I). This material is noteworthy in that both 
of these models fit satisfactorily the increase in the observed number 
of quadrats at two plants per quadrat. As can be seen, the simple 
Poisson series can not possibly fit this particular type of distribution. 


TABLE I 


Distribution of Armeria maritima 


Plants Number of Quadrats 
Per en oe Calculated 
Quadrat Poisson* Neymant | Thomas? 
Ope ener 57 20.6 54.9 57.0 
1 Ateae the 6 3225 7.9 6.0 
VL ts one 12 ray | 10.1 10%1 
S eeteaershs 5 tsyD 9.0 9.5 
Awe.) eaten 5 5.4 6.5 tou.6.5 
WS, fag: Saeee 5 1 BPAY 4.3 4.1 
Greer cege e 7 0.5 2.8 2.6 
t PRUNE 1 (Oe 1.8 1.6 
Secrest is 0 0.0 Dhak 1.0 
OF SN Pe 1 0.0 0.9 0.6 
1 Ossecn wens 1 0.0 0.4 0.3 
11 and over 0 0.0 0.3 0.7 


Total 100 100.0 


* Poisson distribution (Thomas, 1949). 
+t Neyman's Type A contagious distribution (Archibald, 1948). 
* Thomas' double Poisson series, Method II (Thomas, 1949). 


— 
So 
oO 
oO 


100.0 


The number of "plants per cluster" is one of the parameters of both 
these mathematical models. For the Armeria data the results are: 
Neyman, 2.39; Thomas, 2.68. These quantities are of the same mag- 
nitude and can be logically related to a small degree of clumping of the 
plants. In other words, the mathematical model appears to possess 
some relation to reality for this example. For Plantago maritima 
(Archibald, 1948; Thomas, 1949) a similar excellent agreement was 
found between the calculated series and the quadrat frequency data 
The number of "plants per cluster" was found to be 2.16 (Thomas), 


i 
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again a reasonable value. Thomas pointed out, however, that although 
her mathematical set-up was so designed that the parameters of the 
distribution might be capable of physical interpretation, sufficient data 
were not then available to test whether the distribution parameters do, 
in fact, measure the actual number of "plants per cluster" and ''clusters 
per quadrat."' 

A plot of the actual position of every plant of a given species in an 
area usually reveals some type of clumping or clustering which is quite 
evident to the eye but rather difficult to establish rigorously since the 
cluster boundaries are seldom distinct and it is difficult to tell where 
one cluster ends and another starts. An excellent example is shown by 
the position plots of several species in an old-field community on the 
Edwin S. George Reserve of the University of Michigan, Livingston 
County, southeastern Michigan. Cain and Evans (1952) have mapped 
and plotted the distribution patterns of Solidago rigida, Liatris aspera 
and Lespedeza capitata, as they occurred within an area of approximate- 
ly two acres in this community. The number of plants in all the 7640 
one-meter-square quadrats believed to be a valid part of the community 
are studied below. Casual observation of the position plots suggests 
that Liatris is more definitely clumped than Solidago and that Lespedeza 
is still more clumped. 

The object of this study was to determine whether the clusters of the 
contagious series have any relation to these major clumps which are 
evident to the eye. The results strongly suggest that the clusters of 
these series are more closely related to the minor clumps or groups 
of plants that are immediately adjacent rather than to the evident clumps 
(Cain and Evans, 1952), and that one must look elsewhere for a simple 
statistical measure of the degree of major clumping which is evident to 
the casual observer. Such a measure may be related to the average 
spacing parameter described by Dice (1952). 
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CONTAGIOUS DISTRIBUTIONS APPLIED TO OLD-FIELD DATA 


The observed frequency data for the three populations under consid 
eration and the goodness of fit of the double Poisson series of Thomas 
to the values are shown in Tables II to IV. The values have been cal- 
culated by the two methods suggested by Thomas, the moment solution 
(her Method I) and the maximum likelihood solution (her Method II). 

By a mathematical coincidence the second method leads to an exact fit 
of the values for 0 and 1 plant per quadrat. Similar comparisons with 
Neyman's series are shown in Tables II andIV. The calculated values 
of the Neyman series in Table II are more strictly comparable to Meth- 
od I, the moment solution of Thomas, since the constants of the Neyman 
series were also obtained by the method of moments. The application 
of the method of maximum likelihood to Neyman's Type A contagious 
distribution has been described by Shenton (1949). Unfortunately, the 
calculation is involved and requires the usual moment solution as an 
initial approximation. Table IV includes a comparison of the results by 
this method for the Lespedeza population. 

The three species provide examples of an excellent fit (Solidago), a 
fair fit (Liatris), and a very bad fit (Lespedeza), but in all cases the 
size of the clusters of the Thomas model lies between 1 and 3 plants 
(average for the moment solutions, Method I: 1.6; average for the 
maximum likelihood solutions, Method II: 1.3). The size of the clusters 
of the Neyman model are 0.23 (Solidago) and 1.76 (Lespedeza) for the 
moment solutions. These values suggest that the Solidago distribution is 


TABLE I 
Distribution of Solidago rigida 


Plants Number of Quadrats 
Per ie Calculated 
Quadrat Thomas I Thomas II Neyman 
; ail Te Bean eects 7229.0 7229.7 
ace ae 356.0 354.2 
ea my 50.0 50.1 
ree me 5.4 
Sats hacia tc 4 0.5 
5 or over .. 0.0 0.1 
Total .. 7640.0 
Plants per 
clusters || eee 4 j 
Clusters per i ie: 
quadratiie roc . 
Plants per ; Lens 


quadrat |... 
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TABLE II 


Distribution of Liatris aspera 


Plants Number of Quadrats 


ee : Calculated 
uadra 


Plants per cluster. . 
Clusters per quadrat 
Plants per quadrat. . 


TABLE IV 


Distribution of Lespedeza capitata 


Number of Quadrats 


Plants Calculated 
ee Observed Moment Solutions Max. Likelihood Solutions 


Thomas I Neyman Neyman* 


Thomas II 


OR erate, sissy las , 
Dt. a SLCC a en eS 219.6 
Deen RAMEE Gre ich s0'e as 140.8 
Bi Fuck © CO Ree ae 61.6 
AN. ar ae ae aera 210 
mesa ar ots: ie voltat 6! 6.2 
Chit) Gyecr o ROMCae at Oar ORREe iba 
Mitte circ katie cay vio aris tats 0.5 
Caaermey abeuirbriet sre: iol oMs 0.1 
©) So. od @ 618 ahs ORS Ch |i men a al fy AK se maT a A) ee | 
HO eee ae Tee 2 eG fe eel ccccacee ff pw | Basses 
TEL eo abe. citay Sue erieoudioditcdl 6 WC SUNS et (eae came all eal Cel aia ee ee Ty 
DP ROLEONC Tae Aen te eh. Vile seems O) genire s)he 8 OP Tt cathe ee TSF | tee | S(t eaBeeks 
POCA cogtsh Gs eae 7640.0 
Plants-per cluster’... | s.2.... 1.252 
Clusters per quadrat] ....... 0.0853 
Plants per quadrat. .| ....... 0.1068 


*Shenton's method. 
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more regular than random and that the Lespedeza is somewhat clumped 
and point up the difference between the results from the two contagious 
series models. It is evident that these clusters are very much smaller 
than the major clumps which are apparent to the eye on a casual obser- 
vation of the plants themselves. If these species had been randomly 
distributed the number of plants per cluster would have been one, cor- 
responding to the Poisson distribution. The divergence of this parameter 
from unity shows that the plants are not so distributed, but the nature 
and type of nonrandomness is not the same as the expectation suggested 
by the general appearance of the charts. 

It is probable that the minor clusters of the mathematical models are 
only related to very local clumping conditions. This conclusion is partly 
confirmed by an examination of the position plot for Solidago, where a 
large number of very small clumps are evident; the major clumps are 
distinctive and easily differentiated from these local aggregates (see 
also Cain and Evans, 1952). 


MEASURES OF CLUMPING 


A number of authors have considered the general problem of measures 
of clumping without any marked degree of success. Much of the recent 
literature on phytosociological characters has been summarized by Curtis 
and McIntosh (1950), who consider at length the meaning and meanings 
of the fundamental terms such as frequency (proportion of quadrats con- 

- taining species to total quadrats), density (individuals per unit area), 
mean area (area per individual), and abundance (density divided by fre- 
quency). They did not make a detailed study of nonrandomness for their 
synthetic population, although they discovered that the measures of non- 
randomness proposed by McGinnies (1934), Whitford (1949), Blackman 
(1942), and Fracker and Brischle (1944) were all dependent on quadrat 
size when the species were actually nonrandom, 

It must be kept in mind that while disagreement of these measures with 
the expected values from the Poisson distribution certainly indicates some 
kind of nonrandomness of arrangement, agreement does not necessarily 
imply that the distribution of the plants is random since it is simple to 
devise non-Poisson distributions such that the measures are the same 
as those for the Poisson. 

The values of V/X and Ky, Kp, and,Ky in Table V provide these 
four measures of nonrandomness for the one -meter-quadrat data on the 
three old-field species under consideration. These measures (and the 
contagious series results) show that the three species have the following 
types of distribution: 
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Index Most Random to Least Random 
PE OMA S Wats weneucdeaass, Solidago > Liatris > Lespedeza 
Blackman* ...... Solidago > Liatris > Lespedeza 
McGinnies ...... Solidago SS Liatris SS Lespedeza 
ica CKe 1a ert, eens Solidago > Liatris > Lespedeza 
iW Ditton disee .s se. Solidago > Lespedeza > Liatris 


ABS 


The new measure of spacing introduced by Dice (1952) leads to the 
Same conclusion as the measures from the contagious distributions. 

These conclusions can be compared with the general appearance of 
the community as judged by eye from the plots: 


Solidago > Liatris > Lespedeza 


It is evident that for these species most of these measures arrive at 
the same conclusion as the visual estimation. The anomalous result 
from Whitford's index is not surprising in view of its sensitivity to 
quadrat size (Curtis and McIntosh, 1950). There is also no definite 
value for complete randomness. This is in sharp contrast to the other 
measures, 


TABLE V 
Parameters of the Sampled Plant Populations 


Lespedeza 


Mean x or density D....| 0.06191 0.04188 0.10681 
Wiariance.V olin. @ os. . 0.07640 0.07285 0.29435 
Frequency F (per cent) 3.102 6.047 
Abundance A=... 2... : 1.350 1.766 
Expected density* d..| 0.05530 0.03151 0.06238 

Measures of randomness! 

Variance/mean V/x .. 2.756 
Kaas Did) ey eeate tee oe eral 
Kp =(D - d)/d? Sa ae 

Ky = 100D/F2 ...+-- 


* Calculated from the frequency. See Fracker (1944, Table I). 
t+ v/z Blackman (1942), etc. Random: V/x = 

Ky McGinnies (1934). Random: Ky = 

Ky, Fracker (1944). Random: Kr =0 

K, Whitford (1949). Random: Ky is small (0 to 0.3) 
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SUMMARY 


The contagious distributions of Neyman (1939) and Thomas (1949) 
were fitted to quadrat frequency data on Solidago rigida, Liatris aspera, 
and Lespedeza capitata in an old-field community. The goodness of 
fit was excellent, fair, and very bad, respectively. The size of the 
contagious distribution "clusters" of the mathematical model, 1 to 2 
plants per cluster, had little relation to the evident major clumping 
obvious to the eye, but is more likely related to local clumping effects. 

A number of different measures of aggregation or clumping which 
have been suggested in the literature are discussed. None of these 
indicates the actual size of the clumps. This seems to be a new prob- 
lem demanding further research. The relative amount of major clump- 
ing of the three species noted by observation was compared with that 
shown by these different measures and was found to be the same with 
Lespedeza most clumped, Solidago least clumped, and Liatris in between, 

An appendix presents and exemplifies a method for the calculation of 
the individual terms of the contagious series by means of direct cal- 
culation of polynomials instead of by the use of summations. 
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APPENDIX 
Calculation of Terms of Contagious-distribution Series 


It is possible that the ecologists who might have tried Neyman's or 
Thomas' series as a possible mathematical model for quadrat frequency 
data were deterred by the mathematical complexity of these series. 

This section shows that it is possible to evaluate the individual calculated 
terms of these series by the use of simple polynomial functions of var- 
iables which can be determined directly from the experimental data. 
Both series require as starting values the mean and variance of the 
observed distribution. These are obtained by the usual standard pro- 
cedures which are described below. 


Mean and Variance 


Let f be the number of quadrats which contain X plants per quadrat 
(X =0,1,2,....), then the mean or average number of plants per quad- 


rat is: 
My = Z(€X)/ =f 


The summations » extend over all the values of X. 
The corresponding variance or second moment of the distribution may 
be most accurately and usually also most easily obtained from the relation: 


_ 2 (EK?) _ (2 (fx) ) 
Pak Ganeiye bi 
As an example consider the data on the distribution of Solidago rigida 


in the old field (Table VI). 


TABLE VI 


Calculation of Mean and Variance 


Plants Per 
Number of Quadrats Quadrat 


f x {X £X2 
7229 0 0 0 - 
356 1 356 356 
49 2 98 196 
5 3 15 45 
1 4 4 16 
~4 


0 > 0 0 
Total 7640 da 473 613 
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The summations are: 


xf = 7640 Dix = 473 Six°’s=.613 


The mean is: 


m4 = 473/7640 = 0.061911 


The variance is: 


613 / 473)? _ 
Ug = 7640 = or! = 0.076403 


Neyman's Type A Contagious Distribution 

Beall (1939) has shown that calculated values of Neyman's (1939) 
"Type A Contagious Distribution" can be evaluated by a simple method 
involving, for each term, the summation of a linear function of the pre- 
ceding terms of the distribution. This section presents a direct solution 
such that each calculated value of the series can be obtained from a 
polynomial in a new variable without the necessity of calculating the 
preceding terms. The relation was derived by direct substitution in 
Beall's formulas. 


Let m, = Neyman's parameter, proportional to the number of 
"clusters" per quadrat 


m, = Neyman's parameter, proportional to the number of 
plants per "cluster" 


4 = mean number of plants per quadrat 
[ty = variance 


nj = calculated number of quadrats each of which contains 
i plants 


n = total number of quadrats, the sum of the n, 


Neyman has shown that m, and m, may be estimated from the mean 
TA and variance », of the observed values by the relations: 


M2 =(pyo - 14)/ pt 


m, = -Hi/m) = (uf)? / (uy - pt) 
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This is the solution obtained by equating the moments of the theoretical 
distribution and the moments of the observed distribution, The maxi- 
mum likelihood solution has been given by Shenton (1949). This method 
is fairly complicated and requires the moment solution as part of the 
calculation. Shenton's paper includes a worked out example. Beall's 
relations for the calculated values are: 


Pei 
ny = ne m, (1-e ) 
_mime™ EE mt 
ares an i+1 a t!  i-t 


The present relations use a new variable, U, which is defined as follows: 


U=me 2 
The value of n, is: 
ny) =mMm,n U 
The higher terms n,, n,,.... are obtained from expressions such as: 


ny =T oP (1 + U) 


2 
Ns = Ties? (1 + 3U + U?) 


=n, (1+ TU + 6U? + US) 


n, = 
m,' : 
ns =m <7 (1+ 15U + 25U? + 10U3+ U4) 
5 
n, =n, a (1 + 31U + 90U2 + 65U3 + 15U4 + US) 
that is 
if 4 
Se m) : (polynomial in U ) 
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The coefficients of the polynomials with the appropriate values of the 

factorials (i!) are shown in Table VII. The use of these values witha 

table of powers of m, and U permits the direct calculation of any nj - 


Thomas' Double Poisson Series 


The double Poisson series of Thomas (1949) is also presented as a 
summation. 


where m = Thomas' parameter, the number of ''clusters" per quadrat 


X = Thomas' parameter, one less than the number of plants per 
Woelusters 


n, = calculated number of quadrats each of which contains i plants 
n = total number of quadrats, the sum of the n,. 


This equation may be written in the equivalent form which facilitates 
the direct calculation: 


i-r 


Tey! 
-m e r 
pines 1 r! (i-r)! : 
C= 


i-r -A\r 
) 


since the first factor after the summation term is independent of m and 
which must be estimated from the observations. Thomas gives two 
methods of making this estimation. Method I involves equating the 
mean }, and variance of the sample to the mean and variance of the dis- 
tribution. (This is the ''moment" solution similar to that used by Ney- 
man for the other series and described above. ) 


Let g=4H5/p4 
then 2’ =g-34 /g2 -2E95 


and = m= Hy/(1 +n) 
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Method II, the maximum likelihood estimate, involves only the observed 
Nn, and n,: 


tees =n, /n 


e =n, /n,m 


Once m and X are established the calculated values of nj can be 
computed by means of simple polynomials in an analogous manner to 
the treatment for Neyman's series above. A new variable V = me*/, 
is introduced. Direct substitution leads to the following relations: 


-m 
Nh = ne 
-m -2r 
ny = nme e 


Ny =n, 3; (2+V) 


=n, 37 (3 + 6V + V?) 


=) 
| 


3 
n, =n, 4; (4 + 24V + 12V? + V9) 
4 
n, =n, 4 (5 + 80V + 90V? + 20V? + Vv‘) 


n, =n, 2 (6 + 240V + 540V? + 240V? + 30V! + V5) 


that is: 


1 i-] 
mi = —TZr (polynomial in We ) 


The coefficients of the polynomials with the appropriate values of the 
factorials are given in Table VIII. The use of these values with a table 
of powers of A and V permits the direct calculation of any n;, without 
the necessity of considering the previous terms. 
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Example: As an example of the use of these polynomials the cal- 
culation of the terms in the Neyman's Type A series for the Solidago 
rigida data is shown below. 
ui = 0.061911 
LL, = 0.076403 


m, = (0.076403 - 0.061911)/0.061911 
= 0.23408 


m, = 0.061911/0. 23408 = 0. 26449 


e 2 = 0.79130 


m, (1 - e 2) = (0.26449) (0. 20870) = 0.055199 


-m, 
em (l-e “*) _ 9 94630 


n, = (7640) (0.96430) = 7229.73 
U = (0. 26449) (0.79130) = 0.20929 
n, = (0. 23408) (7229.73) (0.20929) = 354.19 


The remaining terms are evaluated using the polynomials. First 
a table of powers of U and m, is prepared: 


Factorial 


U =0.20929 m, = 0, 23408 2 
U2 = 0.043802 m,” = 0.054793 6 
U® = 0.0091674 m,* = 0.012826 24 
Ut = 0.0019186 m,*= 0. 0030023 720 

etc. etc. etc 


It is convenient to jot down the appropriate factorial opposite the prop- 
er power of m, (as has been done above) since it will not only be 
close to the place where it will be used but also serves to locate the 
proper polynomial in Table VIII. 


No. 53 MEASURES OF PLANT AGGREGATION 15 


The calculation of n; serves as an example of the use of the poly- 
nomials: 


n, = (354.19) (0.0030023/120) (1 + (15) (0. 20929) 
+ (25) (0.043802) + (10) (0.0091674) + 0.0019186) 


= 0.047 
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